Energy Flow in Acoustic Black Holes 



K. Choy, T. Kruk, M.E. Carrington and T. Fugleberg 

Department of Physics. Brandon University, 
Brandon, Manitoba, R7A 6A9 Canada 
and Winnipeg Institute for Theoretical Physics, 
Winnipeg, Manitoba, Canada 

J. Zahn, R. Kobes, G. Kunstatter 

Department of Physics, University of Winnipeg, 

Winnipeg, Manitoba, R7A 6A9 Canada 
and Winnipeg Institute for Theoretical Physics, 
Winnipeg, Manitoba, Canada 

o . 

D. Pickering 

Department of Mathematics, Brandon University, 
O ' Brandon, Manitoba, R7A 6A9 Canada 

o" 

We present the results of an analysis of superradiant energy flow due to scalar fields incident on 
an acoustic black hole. In addition to providing independent confirmation of the recent results in 
Q, we determine in detail the profile of energy flow everywhere outside the horizon. We confirm 
explicitly that in a suitable frame the energy flow is inward at the horizon and outward at infinity, 
as expected on physical grounds. 
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The superradiance effect for a rotating black hole describes the situation where fields incident on a black hole have 
a reflection coefficient greater than one, and energy is extracted from the black hole. This theoretical prediction is 
O" 1 extremely hard to test on remote astrophysical objects. In order to study this and other phenomenon of black hole 
physics, it is useful to have a black hole analog which can be studied in the laboratory. The most well known black 
bi). hole analog was developed by Unruh who showed that many aspects of black hole physics could be reproduced in 
supersonic fluid flows [y . The acoustic black hole (or dumb hole) consists of fluid flowing inward toward a fluid sink. 
If the rate of the fluid flow inward becomes greater than the speed of sound in the fluid at some radial distance from 
the sink, sound waves cannot escape from this event horizon. Since Unruh's original proposal, many different kinds of 
analogue black holes have been proposed Q. At present, the Hawking temperatures associated with these analogue 
black holes are too low to be detected, but the situation is likely to change in the near future 0|. 

In this paper we study the rotating acoustic black hole first considered in [3] and recently studied in detail in 
0. The fluid is characterized by an irrotational velocity potential that has a sonic horizon and an ergoregion. The 
corresponding velocity field is given by: 

A B - 

v = r H 4> (1) 

r r 

where A and B are positive numerical constants. We will work in dimensionless variables obtained by setting c = A = 1 
(where c is the speed of sound) . The solutions to the equation of motion for fluctuations in the velocity potential have 
the form: 

■$(t,r,cj)) = R(r)er ^UJt e ^m 'l , (2) 

where to is real and positive and m must be an integer for ^{t, r, 4>) to be single valued. One defines a different radial 
function, H(r): 
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and introduces the tortoise coordinate: 



r* = r + — In 



r - 1 



r + 1 
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which ranges from [—00, 00] as r ranges from [1, 00] (note that r = 1 is the event horizon). We obtain the differential 
equation Q: 



d 2 H(r*) 
dr* 
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H(r*) = 0. 



(5) 



We solve these equations numerically using the asymptotic solutions as boundary conditions. Enforcing the condition 
that the only physical solution at the event horizon is the ingoing one, we obtain: 



H(r*) = Ke lu)r +e~ lujr for r* -> 00 
H{r*)=Te- l{u> - mB)R * for r* -► -00, 

where 7?. and T are the amplitudes of the reflected and transmitted waves and are related by the condition 0, Q : 

mB 



(6) 
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For u> < mB the reflection probability \1Z\ has magnitude greater than 1. Superresonance occurs when solutions 
exist to © that approach the asymptotic forms © for parameters that satisfy this condition. 
The potential function in square brackets in © is shown in Fig. Q for different values of B. 
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FIG. 1: Graph of V(r) vs. r for values of B shown in the legend. 
Our numerical solutions to © agree with those of Results for the reflection coefficient are shown in Figs il'il.'it . 
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FIG. 2: Graphs of \R\ 2 vs. lo for m — 1 and various values of B as shown in the legends. 
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FIG. 3: Graphs of \R\ 2 vs. lo for m = 2 and values of B shown in the legends. 



The superresonance effect is significant even for small values of B, and can be larger than 75% for large values of 
B. For m = 2 the superresonant effect is still significant but decreases with increasing m. Note that superresonance 
does not occur for negative values of m. 



II. SIZE OF THE SUPERRADIANCE EFFECT 



The size of the superresonant effect is somewhat unexpected considering that the superradiance of rotating astro- 
physical black holes is always less than ~ 0.2% It has been pointed out in that another fundamental difference 
between acoustic black holes and astrophysical black holes is that there is, in principle, no upper limit on the rota- 
tional velocity of the acoustic black hole. Related to this is the observation that in the acoustic black hole there is 
no limit on the size of the ergoregion whereas the equatorial radius of the ergoregion in the Kerr black hole cannot 
be any bigger than twice the radius of the event horizon. The superresonance effect, which is related to the Penrose 
process, should have some connection to the size of the ergoregion and, as we will discuss later, the region of negative 
effective potential energy. Neither of these observations, however, account for the difference in the magnitude of the 
superresonant effects. The superresonance effect in acoustic black holes is significant even when the radius of the 
ergoregion is only slightly greater than the event horizon. 

Our results indicate the existence of a maximum value of the superresonance effect as a function of B. (We note 
that the theory behind this analysis is not valid for arbitrarily large values of B. The acoustic metric is derived using 
hydrodynamic approximations, but if the angular component of the velocity is large enough the dispersion relation for 
the fluid will invalidate the assumption of a hydrodynamic model. It is estimated that the theory is valid for B -C 7 
Q). The results are shown in the left panel of Fig. <@}. 




FIG. 4: The left panel is a graph of the maximum value of \R\ 2 versus B for m = 1 using units where A = c = 1. The right 
panel is a graph of |i?| 2 versus u> for m = 2 and 6 = 10 

We find that |7?.| 2 never rises above 2 and is a monotonically increasing function of B. 

It is interesting to note that the existence of a limit to the size of the reflection coefficient can be seen explicitly 
by using a different set of variables. Defining 8 = B/u> in © we find that the reflection coefficient asymptotically 
approaches its maximum value for fixed m and 6 (see right panel of Fig. QJ. These variables were used in the analytical 
analysis of Ref. (which did not consider the limit of large uj). 

Ignoring for the moment the restrictions imposed by hydrodynamics, there is no apparent reason that \T^\inax cannot 
be arbitrarily large, for arbitrarily large values of B. It is interesting to verify that one can reproduce this asymptotic 
behaviour in simple model that can be solved analytically. We consider the general equation 



d 2 H 
dx 2 



V(x)H = (8) 



and use a potential of the form: 



){ x < —a 

V(x) = ^ uj 2 —a < x < a 0J3 = aw\; ui 2 — P^i (9) 
;| x > a 

We use a solution of the form 

H(x) = Aexp(iujix) + B exp(— iuj\x) x < — a 

— C exp(iu> 2 x) + Z?exp(— iuj 2 x) — a < x < a 
= E exp^a^x) x > a 

(10) 

and demand continuity of H(x) and it's derivative at x = ±a to obtain |i?| 2 + a\T\ 2 = 1 and 



(1 - a) 2 cos 2 (2w 2 a) + (j3 - §) 2 sin 2 (2w 2 a) 



\T\ 2 



(1 + a) 2 cos 2 (2w 2 a) + (0 + f ) 2 sin 2 (2a; 2 a) 
4 

(1 + a) 2 cos 2 (2w 2 a) + (f3 + § ) 2 sm 2 (2uj 2 a) 



\R\ 2 + a\T\ 2 = 1 



For (3 ^ 1 superresonnance occurs for a < 0. The analytic results for this simple model clearly indicate the existence 
of a maximum value for |i?| 2 (see Fig. EJ, in analogy with the numerical results for the acoustic black hole. Further 
analysis is required to determine whether or not this is a generic property of superresonance. 
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FIG. 5: Left panel is a graph of \1Z\ 2 versus a with 2au2 = 10. From the top the curves have a — .5/3, a — A/3, a = .3/3 and 
a = .2/3. Right panel is same graph with 2aoj 2 = 5. 



III. ENERGY FLOW 



Our numerical results, obtained concurrently with those of confirm the existence of superresonnant solutions 
for acoustic black holes. If these mathematical solutions correspond to a physical effect, it should be possible to verify 
that there is a net outward flow of energy from the black hole whenever the condition \TZ\ 2 > 1 is satisfied. Moreover, 
for both superresonnant and non-superresonnant solutions, we should find that the energy flow asymptotically close 
to the event horizon is inward, since energy should not be able to escape from within the event horizon. These points 
are well understood in astrophysical black holes j3| but have not been studied in detail in the present context. The 
purpose of this section is to verify that the numerical solutions to the differential equations found in this paper and 
in [5j do indeed give rise (in a suitable frame) to energy flow that is outward in the region asymptotically far from the 
event horizon, and inward in the region asymptotically close to the event horizon, as expected on physical grounds. 

Consider the asymptotic form of the solution very close to the horizon. From © we have: 



H(r*) = T expjicrr*} where a = —{uj — m B), 



(11) 



In the region of parameter space where superresonance occurs (^^p > 1) we have a > and thus, according to an 
observer at infinity, the solution at the event horizon is an outgoing wave with phase velocity given by: 



a . mB 

Vphase \*- )- 



(12) 



Naively, this result might suggest a connection between a positive outgoing phase velocity asymptotically close to the 
event horizon, and the positive energy extracted from a superresonnant acoustic black hole (see Appendix B of 10]). 
A similar result was found in ja| where it was claimed that, for a light pulse with superradiant frequency incident on 
a rotating astrophysical black hole, the energy flowing out of a surface very near the event horizon is positive and 
only slightly less than the energy flowing out of a surface very far from the black hole. However, as mentioned above, 
outward energy flow arbitrarily close to the event horizon essentially corresponds to energy coming out of the event 
horizon. Thus, this naive argument, and the results of seems to suggest unphysical solutions. 

However, as is well known, the direction of the energy flow is governed by the group velocity, and not the phase 
velocity. In this case the group velocity is given by [llj : 



da 

V group — — t 



(13) 



which means that, at the event horizon, energy is flowing inward at the speed of sound, in apparent contradiction to 
statements made in Q (see endnote [lap. 

In order to analyse the energy flow of our solutions, we note that the energy flux across an arc dS of a circle of 
constant r with unit normal vector r" is given by: 



dE = T al3 t a rp dS = T 0r dS. 
where T Q/3 is the stress-energy tensor of the massless scalar field and is given by: 



(14) 



(15) 
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In terms of the numerical solutions for H(r*) the T 0r component of the stress energy tensor is: 



rpOr 



1 — 4m 



2 1 4(Sm— r u) 



H(r*)H*{r*) H '(r*) {H'(r*)Y 



8rYo 2(r2 - l)pg 

i(Bm-r 2 uj) {H(r*) (#>*))* - ff'(r*)IP(r*)} ff(r*) (if'(r*))* + H'{r*)H*(r*) 



2r(r 2 - l)pl 



(16) 



We use the numerical solutions of the differential equation to study the energy flow in both superresonant and 
non-superresonant cases. In non-superresonant cases we find that, as one would expect, the energy flow is inward for 
all values of r. An example is shown in Fig. ©. 
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FIG. 6: Graph of the rate of outward energy flow across a circle of radius r surrounding the acoustic black hole for B = 0.5, 
m — 1 and to = 1. The energy flow is strictly negative indicating that net energy flow is inward everywhere. 

In superresonant cases we find that the energy flow is inward in the region near the horizon and outward far from 
the horizon (see Fig. 0. 
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PIG. 7: Graph of the rate of outward energy flow across a circle of radius r surrounding the acoustic black hole for B — 1, 
m — 1 and U) = 1. Far from the event horizon the rate of outward energy flow is positive which indicates superresonance. The 
energy flow near the horizon (r = 1) is negative indicating that the energy flow is inward. 



As a check of the numerical results, we have confirmed that in the limit r — * 1, the numerical result for the energy 
flow approaches the result that one obtains analytically using the asymptotic solution In both cases the energy 
flow asymptotically approaches a constant value as r — > oo. 

We have also studied how the energy flow, at a fixed value of r far from the black hole, depends on the frequency 
to of the incident wave. The results are shown in Fig. JHJ. 
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FIG. 8: Graph of outward energy flow through a circle at r = 20 as a function of u) for m = 1 and different values of B. As u 
increases, net outward energy flow ceases in each case very close to the cutoff frequency u> — mB, for the superresonance effect. 

Comparing this graph with the rightmost graph in Fig. l|2} we can see that for each value of B, the range of ui for 
which the energy flow is positive is almost exactly equal to the range of w for which [R\ 2 > 1 (which is given by 
u; < mB). We note that the value of u> for which the energy flow changes sign approaches uo = mB, in the limit that 
we look at the energy flow asymptotically far from the acoustic black hole. At finite radius, the frequency at which 
the direction of the energy flow changes is always bigger than lo — mB. 

The radius for which the energy flow becomes positive can be related to a parameter of the well function. From 
Fig. JQ| it can be seen that the well function contains a small negative region. We call the smallest value of r for 
which the potential crosses zero n , and the larger value of r where the potential becomes positive again is called r 2 ■ 
In all cases, the energy flow becomes positive at r > r\. This result is shown in Fig. Jsj. 
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FIG. 9: Graph of the well function for B = 1.6, m = 1 and u = 1 (solid line), and the rate of energy flow as a function of r 
for same parameters (dotted line). The net energy flow becomes outward at a radius greater than the smallest value of r for 
which the potential crosses zero. 

This result is reminiscent of the Penrose process, where energy can be emitted by an astrophysical black hole by 
creation of back to back photons in a region where the effective potential for particles travelling opposite to the 
rotation of the black hole is negative. In the present case, the location of the transition between outward and inward 
energy flow satisfies r > ri, but not r 2 > r > r\ . 

These results confirm that there is a net outward flow of energy from a rotating acoustic black hole in the super- 
resonant regime of parameter space. The energy originates a finite distance outside the event horizon, and the energy 
flow close to the event horizon is inward. We conclude that the numerical superresonant solutions obtained in our 
analysis and that of @ are physical solutions. 
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IV. CONCLUSION 

To summarize, we have obtained, concurrently and in agreement with the authors of 0, a maximum superresonance 
effect that corresponds to an increase in the reflected amplitude with respect to the incident amplitude that is on 
the order of 100% (apart from any limits imposed by hydrodynamics). An analytic solution of a one dimensional 
toy model with similar properties is also shown to exhibit such a maximum. By analysing the energy flow, we have 
explicitly shown that the properties of the energy flow arising in the numerical superresonant solutions are consistent 
with physical expectations. 
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